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Abstract. Given a path of systems of homogeneous polynomial equa- 
tions ft, t & [a, b] and an approximation Xa to a zero Ca of the initial system 
fa, we show how to adaptively choose the step size for a Newton based ho- 
motopy method so that we approximate the lifted path {ft, Ct) in the space of 
{problems, solutions) pairs. The total number of Newton iterations is bounded 
in terms of the length of the lifted path in the condition metric. 



< 



1. Introduction 

Let us denote by 'H(d) the vector space of homogeneous polynomials systems 

/ : C"+i C", 

^ / = (/i, • ■ • , fn) in the variable z = {zq, zi,..., z„), with degree (d) = (di, . . . , d„), 

so that fi has degree di. 

^ Given a path of systems t £ [a,b] ft £ H^d), and a zero of the initial 

I— I system fa, under very general conditions, the path t ^ ft can be lifted to a 

path i — [ft, Ct) in the solution variety 

> t/={(/,C)GH(,) xC"+i : /(C) =0}. 

00 If we make the additional hypothesis that is orthogonal to (t this path is unique. 

Now, given a sufficiently close approximation Xa to the zero Ca of the initial sys- 
tem fa, predictor-corrector methods based on Newton's method may approximate 
the lifted path (/(, by a finite number of pairs {fti,Xi) G 'H(d) x C"+^, < i < fc. 
These algorithms are designed as follows: first the interval [a, h] is discretized by 
a finite number of points a = io < < • • ■ < ^fc = then a sequence {xt) is 
. . constructed recursively by 

• rH Xo = Xa and Xi+i = Nf^ (xi) 

5_j where Nf^ is the projective Newton operator associated with the system ft. The 

complexity of such algorithms is measured by the size k of the subdivision (ti). If 
we make a good choice for (ti) then k is small and, for each i, Xi is an approximate 
zero of /t . associated with Q. . 

The complexity of such algorithms has been related by Shub-Smale in [14J to 
the length I of the path [ft) and to the condition number of the path {ft,Ct)- 
fjL = m.a'x.a<t<b l^ift, Ct)- The condition number measures the size of the first order 
variations of the zero of a polynomial system in terms of the first order variations 
of the system. For (/, C) £ V it is given by 



M/,C) = 11/11 \\{Df iO |c-)"'diag(v/dI||Cr'-') 
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or oo when rank Df{C,) < n. We extend this definition to any pair (/, z) G 
T-L^d) X C"'^^ by the same formula. Shub-Smale in jl4] give the bound 

k < CD^/H)?, 

where C is a constant, and D = maxd^. 

A more precise estimate is given in Shub [13j . The author proves that we can 
choose the step size in predictor-corrector methods so that the number of steps 
sufficient to approximate the Hfted path is bounded in terms of the length of the 
lifted path in the condition metric: 



L = 
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fJ-{ft,Ct)dt. 



We can find such a subdivision (ti), < i < k, with 

k < CD^/^L. 

Its construction is given by 
to = a, and 
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but, in this paper, some universal constants are not estimated, and no constructive 
algorithm is given. 

In the algorithm we present below, we compute ti+i from ti so that at least one 
of the quantities 



dt 





dft 




dt 




dCt 


c 


dt 



dt 



and 



increases of a given fraction of 



D'/^f^ifuXuY 

Moreover, to allow approximate computations, we introduce a tolerance parameter 
£. This algorithm reflects the geometrical structures used in T^. This structure is 
based on a Lipschitz-Riemannian metric defined in the solution variety V by 



where 



(•,-)y,(/,C)MaC)' 
is the Riemannian metric in V inherited from the usual metric on 



P('H(d)) X P(C"+i). This condition metric is studied in more details in Bcltran- 
Dedieu-Malajovich-Shub 2J and f^, Beltran-Shub |6] and [7j, and in Boito-Dedieu 
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Algorithm Homotopy 

Input: (/0*e[a,f,], ^ xq e C"+\ < e < 1/20. 
Output: 

An integer fc > 1, 

A subdivision a — to < ti < . . . < = b, 
A sequence of nonzero points Xi e C"+^, < « < fc. 
Algorithm: 

« ^ 0; io ^ a; 
Repeat 

Find s e [ti, b] so that 



- 5£>>/V(/.,,n) 



In case there is no such s, make s = b. 

Find s' e [ti, 6] so that for all cr € [ti, s'], 

20e2 J, / ^ £ 

^ ' bD^Mf^) - '^'-"^'''^ - 5i?3/2^(/,,,x.) 

where 

0t.,.(a:.) = ||x,|ri||i?^(x,)|;i Ut^x,) - fa{xM. 

In case there is no such s', make s' — b. 
ti+i = min(s,s'); 
Find ^ x^+i € C"+i such that 

(1-3) (iV,,,^^(x.),x.+i) < 

i ^ i + 1; 
Until ti = b. 
Set k ^ i 

This algorithm has the following properties: 

Theorem 1.1. Given < e < 1/20, a homotopy path {ft)te[a,b] ^(d); o,'>T-d 
an initial point 7^ G C"^"'^ satisfying 

£,3/2 g2 



2 r^\./Li7--u//-u\./u,7—w/ - 2' 

then: 

1. a:o is ara approximate zero of fa with associated nonsingular zero Ca, 

2. Lei (/t, Ct)t(E[a,b] ^6 continuous lifting of {ft)t£[a.b] solution variety 
initialized at {fa,Ca)- If the condition length L is finite, then: 

2. a. The algorithm Homotopy with input (e, {ft),Xo) stops after at most 

k = l + 0.65D^/^e-^L 

iterations of the main loop, 
2.b. For each i = 1 . . . k, Xi is an approximate zero of ft^ with associated zero 

Remark 1.1. • This algorithm is robust: it is designed to allow approximate 

computations. 

• For e = 1/20, the computations in \l.l\ and 1.2 have to be exact. 



• The hypothesis "the condition length L is finite" holds generically for 
paths in the solution variety V and, consequently, Theorem holds for 
"almost all" inputs ift)tGla,b]^ ^o- 



4 



JEAN-PIERRE DEDIEU, GREGORIO MALAJOVICH, AND MICHAEL SHUB 



• We reach the same complexity as in Shub |13| . 

• We require the regularity for our path in the space of systems but, in 
fact, an absolute continuous path t € [a, b] ^ ft & ^(d) *s sufficient to define 
the concepts of length and condition length and to prove Theorem The 
proofs given in section are still valid. 

Let us now mention other approaches to the construction of "convenient subdi- 
visions ". 

A practical answer consists in taking t^^i — ti-\-6t for an arbitrary 5t > 0. If x^+i 
fails (resp. succeeds) to be an approximate zero of fti^i then we take 5t/2 (resp. 
25t) instead of 5t\ see Li [TT]. With such an algorithm we may jump from a lifted 
path t — 7> {ft, Q) G to another one t — > {ftX't) ^ ^- Even if, for each i, Xi is an 
approximate sero of /(. , we cannot certify that the sequence (ft-,Xi) approximates 
the path (/t,Ct)- 

In Beltran [T] , the author presents an algorithm to construct a certified approx- 
imation of the lifted path. It requires a C^+^v path in the space of systems, and 
it has an additional multiplicative factor in the number of steps given in [13j . This 
extra factor is unbounded for the class of paths considered here. 

Beltran's algorithm is studied in more detail in Beltran-Leykin T , which contains 
implementations and experimental results. 

Another important problem, which is not considered here, is the choice of both 
the homotopy path (/*), a <t < b, and the initial zero Cq. Classical strategies are 
described in Li [TI], and Sommese-Wampler [18]. 

A conjectured "good choice" (see Shub-Smale [17]) is the system 

(dlz4'-'zi^0, 

/a(zo,zi, . . . ,z„) = y- Cq = (1,0, ...,o), 

I. dn Zq Zji — 0, 

and a linear homotopy connecting this initial system to the target system /;,. See 
[TTj for a precise statement. This conjecture is still unproved. In Shub-Smale [17] an 
adaptive algorithm is given for linear homotopies whose number of steps is bounded 
by the estimate in Shub-Smale [13]. The algorithm we present here is a version of 
that algorithm adapted to the new context of length in the condition metric. 

Beltran-Pardo [5] use a linear homotopy and Beltran's strategy for the choice of 
the subdivision. They get, for a random choice of {fa, Ca), an average running time 
(J{N'^) where N is the size of the input. 

Buergisser-Cucker in |10_ define an explicit algorithm, called ALH, based on the 
linear homotopy and a certain adaptive construction for the subdivision. They 
obtain the complexity 

217D'/^dr{faJb) / li{ftXtfdt 

J a 

which is not as sharp as the estimate based on the condition length given in |13| 
or to our own estimate. Then, we cite the authors, "ALH will serve as the basic 
routine for a number of algorithms computing zeros of polynomial systems in dif- 
ferent contexts. In these contexts both the input system fh and the origin {fa,Ca) of 
the homotopy may be randomly chosen". See this manuscript for a more detailled 
description. 

Our feeling, based on a series of papers on the condition metric: [5], [3], [5], [7], 
[S], is that a good choice for the homotopy path {ft) and the initial zero (a will 
induce a lifted path {ft,Ct) close to the condition geodesic connecting (/a,Ca) to 
{fb, Cb)- We are far from accomplishing this task. 

Our paper is organized as follow. Section [2] recalls the geometric context and 
contains the main definitions. In section [3] we study the variations of the condition 
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number x) when we vary both the system / and the vector x. The main 
difficulty is to estimate universal constants which are already present in many papers 
(Shub-Smale |17j . Shub [13] for example) but which are not given explicitely. Such 
explicit constants are necessary to design an explicit algorithm. Section [4] contains, 
in the same spirit, explicit material about projective alpha-theory. The last section 
is devoted to the proof of Theorem 

2. Context and definitions 

2.1. Definitions. We begin by recalling the context. For every positive integer I € 
N, let Til ^ C[a;o, . . . , a;„], n > 2, be the vector space of homogeneous polynomials of 
degree I. For (d) = {di, . . . , dn) G N", let 'H(d) = Iir^i systems 
/ = (/i ,...,/„) of homogeneous polynomials of respective degrees deg(/i) = di, 1 < 
i < n. So / : C""*"^ — >■ C". We denote by D := maxjdi : 1 < i < n] the maximum 
of the degrees, and we suppose D > 2. 

The solution variety V is the set of points (/, () € T-l(d) x C""*"^ with f{() = 0. 
Since the equations are homogeneous, for all Ai, A2 G C \ {0}, Ai/(A2C) = if and 
only if /(C) = 0. So defines a variety V C P('H(rf)) x P(C"+i) where P(-H(d)) and 
P(C""''^) are the projective spaces corresponding to 'H(d) and C"+^ respectively; V 
and V are smooth. 

Most quantities we consider are defined on ('H(d) \ {0}) x (C"+^ \ {0}) but are 
constant on equivalence classes 

{(Ai/,A2x) : Ai, A2 G C\{0}} 

so are defined on P{7i(^d}) xP(C"+^). This product of projective spaces is the natural 
geometric frame for this study, but our data structure is given by pairs (/, x) G 
(■H(rf) \ {0}) X (C"+^ \ {0}). We speak interchangeably of a pair (/, () in V and its 
projection (/, C) in V. 

Given a path of systems t e [a,b] ft e P('H(d)), and a zero CaP(C"+^) of 
the initial system fa, under very general conditions, the path t ^ ft can be lifted 
to a unique path t — )■ {ft,Ct) in the solution variety V. 

Two important ingredients used in this paper are projective Newton's method 
introduced by Shub in [12 , and the concept of approximate zero. 

For a pair (/, x) G {'H{d) \ {0}) x (C"+^ \ {0}) the projective Newton's operator 
Nf is defined by 

Nfix)^x-Dfix)\-lf{x). 

Here we assume that the restriction of the derivative Df{x) to the subspace orthog- 
onal to x 

x^ = {uG C"+i : (it,x) =0} 

is invertible. It is easy to see that the line throught x is sent by Nf onto the line 
throught Nf{x) so that Nf is in fact defined on P(C"+i). 

Definition 2.1. We say that x is an approximate zero of f with associated zero C, 
(f{C) — 0) provided that the point Xp — Nf(xp^i), xq — x, is defined for all p > 1 
and 

dTiC,Xp)<i-j driCx). 

Here dx denotes the tangential "distance" as in Definition |2.2[ 

A well-studied class of numerical algorithms for solving polynomial systems 
uses homotopy (or path-following) algorithms associated with a predictor-corrector 
scheme. We are given a path (ft), a < t < b, in the space and a root Ca of 

fa- Under certain genericity conditions, the path (ft) may be lifted uniquely to a 
path (/t,Ct) &V,t€ [a,b], starting at the given pair (/a,Ca)- 
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Given an approximate zero Xa of fa associated with Cq, our aim is to build 
an approximation of this path by a sequence of pairs {ft.,Xi) G {H{d) \ {0}) ^ 
^£in+i ^ {0}), 1 < i < fc where, a = < ti < . . . < = h \s a. subdivision of the 
interval [a, 6], and where Xi is an approximate zero of /i. associated with . To 
simplify our notations we let /f. — fi. 

The construction of the suddivision (ti) is given in the "Algorithm Homotopy". 
The construction of the sequence (xi) uses the predictor-corrector scheme based 
on projective Newton's method, studied for the first time in [T3]. This sequence is 
defined recursively by 

Xi+i = Nf^^^{xi). 

In fact, to allow computation errors, we choose Xi^i in a suitable neighborhood of 
^fi+i i^i)- Theorem 1.1 proves that for each i = 1 ■ ■ ■ k, Xi is an approximate zero 
of fti associated with Q. and it gives an estimate for the integer k in terms of the 
maximum degree D, and the condition length of the path (/(, Q), a < t < b. 

Definition 2.2. As it is clear from the context, systems and vectors are supposed 
nonzero. 

(1) *s endowed with the unitarily invariant inner product (see ^) 
where a = (ao, . . . , a„), |a| = ao + • • • + oin, 

/i(^0: ■ • ■ ? ^n) — ^ ^ fija^Q^ ■ ■ ■ X^" , 
|a|=* 

and 

fjii^Oj ■ • ■ ? ^n) — ^ ^ 9i-a^Q "... . 
|a|=<i, 

(2) dR{x,y) is the Riemannian distance mP„(C). 

(3) dp{x,y) = smdii{x,y) = mino^^AeC ^^^jj^jj^^^ is the projective distance. 

(4) dT{x,y) — tandii{x,y) is the tangential "distance''. It does not satisfy the 
triangle inequality. 

(5) One has dp{x,y) < dji{x,y) < dT{x,y). 

(6) When rank Df{x) = n, we denote by 6^ the angle between x and ker Df{x). 
9x — when f{x) = 0. 

(7) 0L,M denotes the angle between the complex lines L and M so that Ox = 

^x.ker Df(x)- 

(8) tp{u) = 1 — 4u + decreases from 1 to Q on the interval [0, (2 — V2)/2]. 

(9) The norm of a linear (resp. multi-linear) operator is always the operator 
norm. 

(10) The condition number 



m(/,x) = 11/11 ||D/(x)|;idiag(v/d;|| 



xW"^-' 



is also denoted fXproj in |14j , and finorm in [S] • 

(11) D = maxdj. We suppose D >2. 

(12) u = D^/^fi{f,x)dn{x,y)/2. 



(13) v^D^/^fi{f,x) 

(14) f3o{f,x)^\\xr^ 



_! 9_ 

l/ll llsll 



where it is assumed that / ^ and g ^ 0. 



\Dfix)\-lf{x)\\. 

(15) 7o(/,x) = \\x\\nmxk>2\\Df{x)\-l^^\\y('^-'l 

(16) ao{f,x)^/3o{f,x)jo{f,x). 

(17) ao = (13 - 3yT7)/4 = .15767 . . . 
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(18) S{f,x) = \\x\\-^Df{x)\^ldi&g{d,)f{x)\\. 

(19) 0^(2^) = \\x\\-^DMx)\-liMx) - fsixM. 

(20) Let a path a < t <h ^ ft £ ^(d) \ {0} ^6 given. We denote by 
the derivative of the path with respect to t, and, for any g € 'H(ci), 



1-911/. 



IWft 



the norm of the projection of g onto the subspace orthogonal to ft divided 
by the norm of ft- The length of (ft) in P('H(rf)) is given by 

rb 



m 



When 1 1 /til = 1 for each t we have 



ft 



dt. 



ft 



m = 



ft 



dt. 



(21) The condition length of the path {ft,Xt) G ('H(d) \ {0}) x (C"+i\{0}), 
a < t < b, is 



ft 



ft 



Xt\ 



1/2 



lJ'{ft,Xt)dt, 



where it is the derivative of the path Xt with respect to t, and where the norm 



ll^*ILt defined as in the previous item with 
for each t we have 



ft 



When \\ft\ 



\xt\ 



1 



L{b) 



b / , 2 

ft 



P^t 



1/2 



H{ft,xt)dt. 



We will also use some invariants related with non homogeneous polyno- 



mial systems. For (d) — (di 



,d„) e W\ let V(d) - nr=i^rf. be the 



set of polynomial systems F — (Fi, . . . , Fn) : C" — > C" in the variables 
X = {Xi, . . . ,Xn), of respective degrees deg{Fi) < di, 1 < i < n. The ho- 
mogeneous counterpart of F is the system f — (fi, . . . , /„) G ^(d) defined 
by 

r / \ di p / -^K 

Ji \Xo , XI, . . . , Xfi ) — '^0^1 1 • • • 1 

\xo xo 

The norm on V[d) is defined by \\F\\ = ||/||. We also let. 

(22) P{F,X)^\\DFiX)-^FiX)\\. 

(23) 7(^,^) =niax,>2||i?F(X)-i^^^||i/(fe-i). 

(24) aiF,X)=PiF,X)^{F,X). 



3. Variation of the condition number 

A necessary ingredient for the proof of Theorem |1.1| is the foUowing theorem 
which gives the variations of the condition number when both the system and the 
point vary: 

Theorem 3.1. Let two nonzero systems f,g G and two nonzero vectors 

x,y £ C"^""^ be given such that rank Df{x)\r^± — n, u < 1/20, and v < 1/20. Then 
rank Dg{y)\y± — n, and 

(1 - 3.805m - v)^i{g, y) < tx{f, x) < (1 + 3.504?/ + v)fi{g, y). 



8 



JEAN-PIERRE DEDIEU, GREGORIO MALAJOVICH, AND MICHAEL SHUB 



Corollary 3.1. LetO < e < 1/4, two nonzero systems f,g G 'H(d)j '"^'^ nonzero 
vectors x,y £ C"^"'^ be given such that u < e/5, v < e/5, and rank Df{x)\x± — n. 
One has rank Dg(y)\y± — n, and 

(1 - e)^i{g, y) < fj.{f, a;) < (1 + e)n{g, y). 

The proof of these results is obtained from the following series of lemmas. 

Lemma 3.1. Let f € ^{d) (^nd x G C"+"'^. We have 

(1) For any i = 1 . . .n, ||/i(x)|| < ||/i|| ||a;||''% so that ||/(a;)|| < 1 when f and x 
are normalized. 

(2) For anyi = l...n, \\Df,{x)\\ < di \\f,\\ so that \\Df{x)\\ < D when 
f and X are normalized. 

(3) Jo{f,x)<^^l{f,x). 

Proof. These inequalities come from Proposition 1, and Theorem 2, p. 267, in 
i- □ 

Lemma 3.2. 1 < y/n = minf x x). 

Proof. Let given a matrix A E m > n, and A = UY,V* a singular value 

decomposition with 

E = diag(ai(A) > . . . > an{A)) e C"^", 

o'i{A) > 0, U E U„, V e Vjn unitary matrices. We define 

k{A) = ar,{A)-^\\A\\F = <Jr,{ArWa,iA)^ + ... + <Jn{Ay 

when cr„(A) > 0, and oo otherwise. We see easily that 

min k(A) = ^/n, 
A 

and this minimum is obtained when ai(A) = 1, 1 < i < n. 
For the case of polynomial systems we have 

= \\f\\\\(diag{d;'/^)Dfix)\,^y'\\ = (diag(d-^/')i?/(x)U.) "\ 

Using the unitary invariance of the norm in that is ||/ o U\\ = ||/|| for any 

U € U„+i, considering a U such that Ueo — x with eg = (1, 0, ... , 0)"^ e C"+"'^, we 
see that 

||diagK-^/')i?/(x)U.|l^<||/||. 

Thus, our previous estimate shows that ^{f, x) > yfn. 

To prove the equality \/n — minj^ a;) we use the unitary invariance of the 
condition number 

2;) = ti{! o U,U*x) 

for any U £ U„+i, and the equality y/n — cq) when fi{z) = y/diZQ^~^ Zi, 
1 < i < n, and eo = (1, 0, . . . , 0)^ G C"+i. □ 

Lemma 3.3. For any f and x one has dji{x, y) < u/\/2n < u/2, and du{x, y)6{f, x) < 
u. 

Proof. We suppose that both system / and vectors x, y are normalized. We 
have by Lemma |3.2| 

253/2 23/2 
u = -^fi{f,x)dR{x,y) > -^VndR{x,y). 
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For the second inequality we have 

dR{x,y)6{f,x) < d«(x,2;)p/(a;)|;idiag(dl/')|||idiag(dy')||||/(a;) 

< dR{x,v)p^{f,x)D^'^\\f{x)\\ 

< dRix,y)fi{f,x)D^/y2 
= u 



thanks to the inequahties 2 < D, \\fi(x)\\ < (Lemma 3.1), and the hy- 

pothesis 11/11 ~ \\x\\ = 1 . □ 



Lemma 3.4. When rank Df{x)\r^± — n we have 

\\Dfix)\-lDf{x)\y^\< 1 + dp{x,y)taii9,. 

Proof. Take u G y-^ and define v — Df{x)\'^lDf{x)\y±u so that 

V = {u + kcr Df{x)) n x^. 

V € a;"'" is the projection of u € y*" along ker Df{x). Let us denote by w the 
orthogonal projection of u onto x-^. See Figure 1. We have = cos9u,w and 




Figure 1. 



\w — v\\ — \\w — u\\ tan^a; — \\u\\ sinOu.w tan^j;, 



so that 

\\v\\ < llulKcos^u^^ + sin 6'„,u, tan 6*2:) < ||u||(l + sin^x.jf tan^x) 
because 0x-i-,y^ = O^^y = max ||^|| ^ ^ 9u,w □ 

u E y-^ 



Lemma 3.5. Assume that \\x\\ — \\y\\ — 1. When rank Df{x)\r^± — n and u < 1 
one has 

\\Df{x)\-lDf{y)\y^\\< ^^-^ +dpix,y)t^n6.,. 
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Proof. Since Df{y) = Df{x) + Y.k>i ^^^^TT^iv " we get 

r)k+i f(\ 

Df{x)\-lDf{y)\y. = Df{x)\7jDf{x)\y.+J2{k+l)Df{x)\-l ^^y' (y-^)%- 



Then, we use \\Df{x)\^l ^(fc+ij^^ II - loif,x)'', Lemma 



3.4 



and Lemma 



3.1 



to obtain 



\\Df{x)\-lDfiy)\y^\\ < 1 + dpix,y) tanS, + ^(fc + l)7o(/, x)''\\x - yf 



< 



(i-uy 



k>l 

dp{x,y) ta,n9x 



Lemma 3.6. When rank Df{x)\j.± — n then 

Une^^5{f,x)<D^/^ti{f,x). 

Proof. We assume that x and / are normahzed. Let y — Df{x)\~lDf{x)x 
be the projection of x onto x^ along ker Df{x) so that ||y|| = tan 0a;. By Euler's 
identity, one also has 

y = Df{xXldiagid,)f{x) 

thus, 

tan 6*3; = S{f,x) 

< \\Df{x)\~ldiag{dy')\\\\die.gidy')f{x)\\ 



Using Lemma 3.1 we obtain 



tan6'a; < D^^^fi{f,x). 



Lemma 3.7. When rank Df(x)\^± = n one has 

\\Dfix)\;lDfix)\y^\\< 1 + Sif,x)dpix,y). 
Moreover, when \\x\\ = \\y\\ and u < 1, we have 



\Df{x)\-lDf{y)U\\< 



1 



+ S{f,x)dp{x,y). 



Proof. The first assertion comes from lemmas [3^ and [3761 The second assertion 



is a consequence of lemmas 3.5 and 3.6 



□ 



Lemma 3.8. Assume that \\x\\ = \\y\\ = 1. When rank Df(x)\^± — n, and u < 
(2 - %/2)/2, then rank Df{y)\^± = n ana 



\\Df{y)\-lDf{x)l.\\< 



Moreover, if u < 1/19, 



\\Df{y)\-lDfix)U.\\ < {l + dpix,y)Sif,y))^-^j^ < 1 + 3.805u. 



Proof. We have 



k>l 



SO that 



Dfix)\-liDfiy) - Dfix)) J2 Dfi^rj^^^^iy ^ x)' 



k>l 
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and, like in the proof of Lemma |3.5 

{1-u) 



Df{xXl{Df{y) - Df{x))\\ < Y^ik + = < 1 

fc>i 



(because u < (2 — \/2)/2) so that, by Neumann's Perturbation Theorem, 

4^ + Df{x)\-l{Df{y) - Df{x))\,^ = DJ{x)\-lD}{y)\,^ 
is invertible (i.e. Df{ii)\^± is invertible) , and 

1 (1-w)^ 



\Df{y)\-}Df{x)l^\\ < 



1 - ffE$ V'(^) 



We will prove the second statement in two steps. First, prove it under the 
assumption that Df{y)\y± is invertible. Then, we remove this assumption. 

The first step goes as follows. Combining the first statement with Lemma [3^ we 
obtain: 

\\Dfiy)\-lDfix)U\\ < \\Dfiy)\-lDf{y)U.\\\\Df{y)\-lDf{x)U\\ 



< {1 + dp{x,y)taTi9y)- 
= (l + dp(x,y)<5(/,2/))' 



■0(u) 

A bound for 5{f, y) is 

<5(/,y) = \Df{y)-lAi^g{d,)f{y) 

< \\Dfiy)\-lDf{x)l^\\\\ \\Df{x)-ldmg{d,)fiy)\\ 

< \\Dfiy)\-lDf{x)l^Mf,xWD 

using |/(?;)| < 1 since ||/|| = 1 and ||y|| = 1. 

Combining both inequations and setting M = \\Df{y)\~lDf{x)\^±\\, we obtain: 

(1 - u]'^ 
M < (1 + Mu) \ , / 

that simplifies to 

^ = . }''f. 3 < 1 + 3.805.. 

The last bound follows from the fact that the numerator and the denominator have 
alternating signs, so the Taylor expansion at zero of the fraction has terms of the 
same sign (positive). Hence, 



1 — 5'U-t-4u^ — u-* 



1 



u 

is an increasing function. In particular, for u = 1/19, this is smaller than 3.805. 

Now, we must prove that Df{jj)\y± is invertible. Let {xt)t£[o,dii(x,y)] denote a 
minimizing geodesic (arc of great circle) between x and y. 

Let W be the subset of all t G [0, ^^(a::, y)] so that Df{xt)\xt is invertible. It is 
an open set, and Q €W. 

We claim that is a closed set. Indeed, let s € W. Then there is a sequence of 
ti £W with ti — >■ s. We know from the second statement (restricted) that 

\\DfixtX'Dfix),^\\ < 1 + 3.805U. 



12 



JEAN-PIERRE DEDIEU, GREGORIO MALAJOVICH, AND MICHAEL SHUB 



Hence, for t = ti, 

h{T) = \\Df{x,)\-lDf{x),^\\l < n(l + 3.805u)2 

The function /i(t) is a rational function of a real parameter r, so its domain is an 
open set and contains s. By continuity, h{s) < n(l + 3.805u)^. Thus, Df{xs)x^ is 
invertible, and s S W. As 14^ is a non-empty open and closed subset of an interval, 
W — [0,dji{x,y)] and Df{y)\ ± must be invertible. □ 



x) < fi{f, y) ( — ^ + 8{f, x)dp{x, y) 



Lemma 3.9. Suppose that rank DJ{x)\x-i- — n, u < \ and y) is finite. Then 

1 

Proof. Suppose that ||a;|| = ||y|| = ||/||= 1. We can bound 

= \\Dfix)\-lDf{y)\y^Dfiy)\-ldmgi^,)\\ 
and we conclude with Lemma 13/71 □ 



Lemma 3.10. When rank Df{x)\x± — n, and u < 1/19 we have 

Kf,y) < (l + 3.805u)/i(/,x). 
Proof. Suppose that ||2^|| = ||y|| = ||/||= 1. We have 

M/,y) = \\Dfiy)\-lDfix)U^Dfix)\-ldiiigi^/d,)\\ < (l + 3.805M)^(/,a:) 
by Lemma [XS} 



Lemma 3.11. Assume that \\f\\ = \\g\\ = 1. Suppose that rank Df{x)\x± — n, and 
V < 1. Then rank Dg{x)\x± = n, and 

\\Dfix)\-lDgix)U\\ <l + v, 
\\Dg{x)\-lDfix)U\\ < 

1 — f 

(1 - v)fi{g, x) < ^{f, a;) < (1 + v)^i{g, x). 

Proof. Suppose that |la;||= 1. One has Df{x)[^lDg{x)\^± — — {I^± — idem) 
and 

4^ - idem = Df{x)\-ldiag{dy^)di^g{d-^/^)D{f ~ g){x)\,^ 
which norm is bounded by (using Lemma 3.1) 

fi{f,x)D'/^f^g\\<v<l. 

This proves the first inequality. Thus D f (x)\~lDg(x)\x± is invertible and the norm 
of its inverse is bounded by 1/(1 — v) (Neumann's Perturbation Theorem). This 
gives 

f,{g,x) < \\Dg{x)\-lDf{x)U.\\\\Df{x)\~ldi^g{{dy'^ < f^^. 
The last inequality is obtained via 

fi{f,x) < \\Df{x)\-lDg{x)U^\\\\Dgix)\-ldiagiid]/'^ < (l + «)Mg,x). 

□ 

Lemma 3.12. Let Ug — D^^^ ^{g, x)dii{x,y)/2. Suppose that rank Df{x)\r^± = n, 
and V < 1. Then 

(1 — v)Ug < M < (1 + v)Ug, 

and 

(1 - v)5{g, x)-v< 5{J, x) < (1 + v)8{g, x) + v. 
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Proof. Suppose that 1. The first double inequaUty is a consequence of 

Lemma [3. 11 1 For the second one, one has: 

= pg(x)|;ii?/(x)U^I?/(a:)|;idiag(d,)/(a;)+i^g(a;)|;idiag(d,)(.g(a:)-/(.T))|l 
By Lemma ^Al\ \Dg{x)\^lDf{x)\^±\\ < 1/(1 --y) so that 

6ig.x)<^-^+f,{g,x)D'/^f-g\\. 
1 — V 



Again by Lemma 3.11 ii{g,x) < fi{f,x)/{l — v) so that 

,^ Sif,x) , tiif,x)D'/^f-g\\ ^ Sif,x)+v 
o[g,x) < \ < — . 

1 — V 1 — V 1 — V 

Similarly S{f, x) — 

\\Df{x)\^lDg{x)UDgix)\-ldiag{d,)gix) + Df{x)\-^diag{d,){f{x) - g{x) 
\\Df{x)\-lDg{x)\,^\\S{g,x) + v. 
. shows that \\Df{x)\'^lDg{x)\^±\\ < 1 + w so that 
S{f,x) < il + v)Sig,x) + v. 

a 



Lemma 3.13. When rank Df{x)\^± — n, u,v < 1/20, then 



< 



Lemma 




1 + 3.805. 



Proof. As before, let Ug = D^^^ iJ.{g,x)dR{x,y)/2. Lemma 
bound Ug < u/{l — v) < 1/19. So we may apply Lemma 3.10 to g instead of / so 
that 

f^{9,y) < {l + 5-805ug)n{g x) 
Then, we bound fi{g,x) by fi(f,x)/{l — v) (Lemma 3.111 and obtain the first in- 
equality. In particular, fi{g, y) is finite. 

To prove the second one we apply Lemma 3.11 and Lemma |3.9| to obtain 

2 + ^{9:X)dp{x,y) \ fJ,{g,y). 



/i(/,2:) < {l + v)^i{g,x) <il + v) 
By Lemma |3.3| and Lemma |3.12| we have 



(1 - %) 



S{g,x)dp(x,y) < Ug < 



1-v 



and we are done. 



□ 



3.1. Proof of Theorem |3.1[ To prove the inequalities (1 — 3.805 — v)fi{g,y) < 
/^(/j x) < {1 + 3.5047i + v)ii{g, y) we use Lemma 3.13| which gives 

B{u, v)fi.{g, y) < fi{f, x) < A{u, v)fi{g, y) 

with 



A{u,v) = {l + v) 



1 



1 - 



u 

2 1 — V 



, B{u,v) = 



l + 3.805y^ 



A{u, w) = l + 3u + w + 'u 



+ 9™^ - 12u^ + Au^v - I2uv + Qv - 2u^ + iu 



(1 — w) (1 — M — w) 

and the last parenthesis is less than 0.504 when u,v < 1/20. 

The function B{u, f ) — (1 — 3.805u — v) is increasing in u and v, and vanishes at 
the origin. □ 
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3.2. Proof of Corollary 3.1, Since u, v < e/5 and e < 1/4 we get u,v < 1/20. 



Thus we can apply Theorem |3 . 1 1 which gives 

i^~s)^i{g,y) < il-3M5u-v)ti{g,y) < ^i{f,x) < {^+3.50'iu+v)^lig,y) < {l+s)ti{g,y). 
□ 

4. Alpha theory in projective spaces 



Theorem 4.1. Let < a < uq = (13 - 3V17)/4 = 0.15767. . . Let f e U^d) and 
X G C"^"'^ both nonzero. If 

-^Kfi^)l3o{f,x) < a, 
then there is a zero C G C'^"*^ of f satisfying: dxix,^) < a{a)l3o(f,x) with 

1 , l-V(l+a)2-8a 

(jia) = — . 

^ ^ 4 4a 

Furthermore, if y — Nf{x), then (i_R(2/, C) ^ ("'(o^) ~ 2;). Moreover, when 

then X is an approximate zero of f corresponding to C, and so does y. 

Proof. The proof below follows the lines of [T7]. We suppose that / and x 
are normalized (||/|| — \\x\\ = 1). We consider the non-homogeneous polynomial 
system, defined for a variable X G by 

F{X) = f{x + X). 

Let us denote by Np the usual Newton operator: 

Nf{X) = X - DF{Xy^F{X). 

Then DF{X) = Df{x + X)\^±. In particular, DF{Q) = Df{x)\^± and we have, 

y^Nf{x) ^ A(j; + iVF(0)),AG C\{0} 

Po{f,x) = /3(F,0) 

lo{f,x) > ^iF,0) 

aoif,x) > a{F,0). 

Since, by Lemma 3^ aa{f,x) < /3o{f, x)fj,{f, x)D^^'^ /2, by Theorem 1, p. 462, in 
[14| . is an approximate zero of F and hence (Definition 1 ibid.) the sequence 
{Xk)k>o defined recursively by Xk+i — Np{Xk), Xo = 0, converges quadratically 
to a zero Z of F. Namely, 

<2-2'=+i||Xi-Xo|l. 

Moreover, by the same theorem, 

II < 1 + o(f , 0) - v/(l + a{F, 0))2 - 8^(j^ 
" 47(F,0) 
Thus, for C = (x + Z) / ||a; + Z||, we can bound 

dTiC.x) < \\Z - XoW < a{aoif, x))l3{F,0) < a{a)Mf,x). 
Again by Theorem 1, p.462, of [H], 



1 - 3a(f , 0) - V(l + a(f , 0))2 - SajF, 0) 

" - Mm 

which implies that 

dniCy) < {a{a)-l)Mf,x). 
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Let us prove that x is an approximate zero. This wiU follow directly from [5], 
Chap. 14, Theorem 1. In order to apply this Theorem, we have to check its 
hypothesis 

dT(C,a;)7o(/,C)<^^. 

3- \/7 



Using Lemma |3.1| this is obtained from 

^3/2 



(4.4) 



rfT(C,^)^M/,C)< 2 



We notice that u = ^^^dii{x,()fJ-{f,x) < ut — dr {x X) f^if , x) < aa{a) < 
0.049f7(0.049) = 0.0518 • • • < 1/19. According to Lemma [STO) 

M/,C)<l-2M(/,a;). 

Hence, we infer (4.4) from: 

driC, x)—^i{f, < 1.2«T < 0.1771 • • • = 
A similar argument holds to prove that y is an approximate root: 

dTiC,y)^KfX) < l.2dTiCy)—fi{f,x) < 1.2a(a(a) - a) < 
□ 

In the following proposition we relate the invariant /3o(/, a;) for an approximate 
zero X to the distance from its associated zero 

Proposition 4.1. Let f E be fixed and x E C"+^ be given. If 

^/3o{f,x)fi{f,x)<a<0m9, 

then, l3o{f,x) < 1.128dT{x,(^), where C is the zero of f associated to x, given by 
Theorem \4-l\ 



Proof. We suppose that both /, x, and C, are normalized. From Theorem |4.1 
drixX) < PQ{f,x)a{a). Hence, 

£,3/2 

u<UT^ -—n{f,x)dT{x,Q < aa{a) < 0.0518. 



From Theorem 3.1 we conclude for later use that 

£,3/2 

UT,c = -^t^if> Odrix, C) < ^^(1 + 3.8057/) 

Now we can bound: 

/3o(/,x) = \\Df{x)\-lf{x)\\ 

< \\Dfix)\-lDf{C)cA\\\Df{0\-lfix) 



< 



+ u]\\Df{C)\-lf{x 



(1 -m)2 

using Lemmas |3.5| and |3.6[ We further bound, as usual 
\\Df{C)\-lfix)\\<\\x~C\\+Y.l^i^f(0\c'D'fiOix^O'l^ 



k>2 



Putting all together, 
/3o(/,x) < 



1 



-dT{xX) < l.l28dT{x,Cy, 



(l-?i)2 1 - u(l + 3.805w) 
the last step is obtained numerically, using u < a<T{a) < 0.0518. □ 
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Proposition 4.2. Assume that 

^3/2 



Mf,x)^lif,x) <a<l/20 



2 

Lety = Nf{x). Then, 

Proof. We assume ||a;|| — 1. Let F : X e x-^ t-^ f{x + X) be the afHne 
polynomial system associated with /. Then, /3(F, 0) = /3o{f,x). Moreover, we can 
scale y = x + Y, ioi Y = Nf{0). By Proposition 3, p. 478 in [11], 

(4.5) l3{F,Y)<^^j^f3o{f,x). 

via) 

Moreover, 

(4.6) Mf,y) < \\Df{y)\~lDfiy)\,^\\\\y\\~'mY). 



Clearly, \\y\\ > 1. It remains to bound the norm of the first term in the rhs of (4.6 1. 
By hypothesis, 

£,3/2 

u = -—fi{f, x)dR{f, x)<a< 1/20 



so Theorem 3.1 implies that 



1 — 3.805a 



In particular, Df{y)\yj_ has full rank, and we can apply Lemma 3.7 and then 
Lemma \3M to bound 

\\Df{y)\-lDf{y)U4 < 1+S{f,y)dp{x,y) < l+D'/'fi{f,y)/3o{f,x) < 1+ " 



Combining with ( |4.5[ ) and (4.61, we obtain: 

My) < ^ (i + ^-i^ ) ML^y ^ i-23aA(/,x) 



5. The homotopy 

The objective of this section is to prove Theorem |1.1| Through this section the 
considered systems and zeros are normalized: /(C) = with ||/|| = \\(\\ = 1. 
Let t e [a,b], and xt £ C"+^ be given with ||xt|| — 1. We suppose that 

(5.7) rank Dft{xt)\^j.± = n, 
and that 

£,3/2 

(5.8) —l3o{ft,xt)fi{ft,xt)<a. 



According to Theorem 4.1 for a small enough, xt is an approximate zero of /(. We 
call (t the associated zero and extend it continuously for s e [t, t'] so that fsiCs) — 0. 
The main difficulty to prove Theorem |1.1| is to transfer the properties (5.7) and 



(5.8) supposed to be true at t — ti onto a similar property at t' — ti+i. Moreover, 
we must show that if Xt is an approximate zero associated to then the same is 
true for f, for a continuous path (.s- For this purpose we study this transfer in a 
general context. 

Through this section, e < 1/6. Let t' > t he given and assume that 

(5.9) m-iit)< ' 

bD^l^lx{ft,Xt) 
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and 
(5.10) 



max 6tA^) < -. 

sG[t,t'] 5D-y^n{ft,xt) 



For any s let us define Xs = Nf^{xt). Notice that Xs is not necessarily 

normalized. 

Lemma 5.1. Let e < 1/6 and set a = Under the hypotheses above, for any 

s €]t, t'], one has 

(1) n{fs,xt) < jhifJ.{ft,Xt). 

(3) ^^-^/3o{fs,Xt)fJ.{fs,Xt) < 0.049. In particular, Xt and Xs are approximate 
zeros of fs associated with (^^ , 

(4) ii{fs,Xs) < J^,^^{ft,xt). 

(5) {1 - e)n{fsXs) < t^ift,xt) < (1 + e)^(/s,Cs)- In particular, Cs is non- 
degenerate zero of fs, and hence s i— )■ Cs in continuous for s G [t,t']. 

(6) Ms.^s) < 1.23aifs,Xt)Mfs,Xt). 

(7) Hypothesis (5.7) and a strong version of (5.8) hold ats: rank Dfs{xs)\x^ = 
n, and 

£,3/2 

-—/3o{fs,Xs)fi{fs,Xs) < 0.128a 



Proof. 1. From equation (5.9), 



\\ft~fs\\<i{s)-i{t)<m~i{t) < 



5D^/^fi{ft,xt)' 



sothat V = VD^{ft,Xt)\\ft~fs\\ < e/5 < 1/20, and Corollary |3.1 [gives (1 - e) Xf ) < 
M(/i,a;t). 



Mfs,xt) = \\Dfs{xtXlfsixt)\\ < 

\\Dfs{xt)\-lDf,{xt)\,^\\ {\\Dft{x,)\-l ifsixt) - ftixt)) II + \\Dft{x,)\-lft{xt)\ 
By Lemma 3.11 ( |5.8[ ) and ( 5.10[ ) we obtain: 



l-v 



1 



2a 



e/5 + 2a 



l-e/5 V5i?3/V(A,x,) Dy^^iift,xt)J {l-e/b)D^l^n{ft,Xty 
For the lower bound, 

Po{fs,Xt) = \\Dfs{xt)\-lfs{xt)\\ > 

-1 ~^ 

(Dfs{xt)\-lDft{xt)l^) {\\Df,{xt)\-l ifsixt) - ftixt)) II - \\Dft{xt)\;lftixt)\\ 

By Lemma [3lT| [5?8l and [SlOl we obtain: 

1 



Poifs,Xt) > 



1 + v 



kA^t)- Mft,xt)) > 



> 



1 



4>t,s{xt) - 



2a 



l + e/5 

3. Combining the the two preceding items, 
(5.11) 



DV^fi(ft,Xt] 



g/io + Q, 

Mfs, xt)fi{fs,xt) < ^ < 0.037931 • • • < 0.049. 



2 '-w.,-^wr-w.,-^w - (l_g)(l_g/5) 

Thus, by Theorem |4.1[ Xt and Xg are approximate zeros of fs associated with Qs 
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4. Using item 2, 

dnixt^Xs) < /30(/s,Xf) < /t( f 

(1 - e/5)D-^/^^i{ft,Xt) 

Thus, 

u = D^'''dR{xuXs)ix{fuXt)/2 < ^, " < e/5. 

(1 - e/5) 

Then, we can use Corollary |3 . 1 1 again to bound 



5. From Theorem |4.H 

dRixtXs) < drixtXs) < Tiaifs,Xt))Mfs,Xt) < 1.0429 ••• /3o(/s, x*). 

Thus, 

£,3/2 

u = —tiift,xt)dR{xtXs) < 0.1978026 • • • e < e/5. 

We bounded in item 1 the quantity v = D^/^ ^{ft, Xt)\\ft — fs\\ < e/5. Hence, by 
Corollary |3.1| again: 

(1 - e)^i{f,Xs) < KftXt) < (1 + e)^l{fs.Xs). 

6. From item 3 and Proposition |4.2[ 



Poifs,Xs) < l.23ao{fs,xt)l3o{fs,xt). 



7. Because fi{fs,Xs) is finite, Dfs{xs) has full rank. From items 1, 6 and (5.11), 



l3o{fs,Xs)fi{fs,x,) < 1.23 



e/lO + a 
{l-eKl-e/5) 



< 0.128a. 



□ 



Recall that our algorithm allows for an approximate computation of the Newton 
iteration. The robustness Lemma below shows that if a point x satisfies (5.8) and 
conclusion 7 of Lemma 5.1 then an approximation y oi x satisfies (5.8) and (5.9). 

Lemma 5.2. Assume that ||/|| = 1 and \\x\\ — \\y\\ — 1. Let a < 1/72 and c < 0.8. 
Suppose that Df(x)\.j.± has rank n, and 



(5.12) 
(5.13) 



^3/2 



M.f,x)^l{f,x)<0.U8a 



dRix,y)fi{f,x) < 



Then, Df(y)^y± has rank n, and 

£,3/2 

(5.14) — 



Dpi{f,x) 



and furthermore, x and y are approximate zeros associated to the same exact zero 



Proof. By using D^/^^{f, a;) > 4 (see Lemma 3.2 and the hypothesis £> > 2) we 
obtain that u < 0.0055 • • • < 1/19. Therefore, Lemma 3.10 implies that 

m(/,2/) < (l + 3.805)M(/,a:) 

and in particular, Df{y)\y± has rank n. 

To estimate Po{f, y), we decompose 

< 



/3o(/,2/)= Dfiy)-^ifiy) 



Df(y)-;^Df{x) 



Df{x)-^\f{y) 
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The first term is bounded by Lemma [3^ 



< 1 + 3.805m. 



Taylor's exansion gives Df{x). \f{y) = 

Dfix)-^\fix) + Df{x)-^^DJ{x){y ~ x) + ^^D f {x)-^\D^ J {x){y - 

k>2 

Taking norms, 

Df{xV\f{y) < /3o(/, x) + S{f, x)\\y x\\ + H^^ - ^ll'7o(/, x) 
' 1 - ||y a;||7o(/,x) 



By Lemma 3.1 :,\\y — a;||7o(/, x) < u. Hence, 

(io{f,x) < (1 + 3.805U) (f3o{f,x) + 5{f,x)\\y-x 



\\y - x\\u 

l-u 



Using Lemma 3.6 S{f,x) < \/D^{f,x). Thus, 

y)lJi{f: y) < (1 + 3.805u)2 f 0.128a + VDfi{.f, x)u + 



l-u 



< (l + 3.805u)2 ( 0.128 + c 

< 0.97q! < a 



cu 



Dfi{f,x){l-u) 



a 



Since a < 1/72 < 0.049, Theorem 4.1 imphes that both x and y are approximate 
zeros of /. As this is also the case for all the points in the shortest arc of circle 
between x and y, the associated zero must be the same. □ 



Lemma 5.3. Assume the Hypotheses of Lemma \5 . 1\ Let < ^ < 1. If furthermore 

l{t') - l{t) > 



5VDiJ,{ft,Xt) 



the 



L{t') - L{t) > 



Proof. 



m-m 



> 



> 



> 



> 



i: 


ifsXs) 


KfsXs) ds 


i: 


fs 




ds 


Kft,xt) /■* 
l + e Jt 


Is 


ds using Lemma [5T|5) 


Kfuxt) 
l+e ' 


-lit)) 












5(1- 










€e 











5£'3/2 ■ 
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Lemma 5.4. Assume the Hypotheses of Lemma \5.1\ and choose e and ^ so that 
20e < ^ < I. If furthermore 



then 



L{t') - L{t) > 



5D^/^fiift,Xt)- 



13i:»3/2 ■ 



Proof. 



UsXs) Kfs,Cs) ds 



> 



> 



> 



Cs 



1 + e 
1 + e 



KfsjCs) ds 



Cs ds using Lemma |5.1[ 5). 



dRiCtXf). 



By the triangle inequality, 

dRiCtXt') > dB{Ct',xt) - dnixtXt) 

We know from Theorem 4.1 that dji{xt,Ct) < (^{ct)Po{ft,Xt)- Lemma [5.1[ 3) says 
that ctQ{ft',Xt') < 0.049 and hence, from Proposition 4.1 we obtain: 

dTixtXt')>Mft',Xt)/1.128. 

Thus, dxixtXt') > {i+e/5)Dml^{ft,xt) T328 Lemma 5.1 [2). We use now the bound 
20e < ^ and the fact that a = to deduce that 



dxixtXt') > w = 



3Ce/20 



1.128(1 + e/5)i?3/2/i(/t,^t) 
Since e < ^/20 < 1/20, we can bound lo < 0.001645 • • • . Of course, dfl(a;t,Ct') > 
arctancj. We may bound arctan(w) > cj arctan' (0.001645 ••• ) = Y+cTooIeiis^'^ ^ 
0.999w. Now, 



L{t') - L{t) > 



tJ-ift,Xt) 



> 0.0775 



0.999 



> 



3$e/20 



2(T(a)a 



153/2 - 131)3/2 



1.128(1 + e/5)D^/^n{f„Xt) DV^^l{ft, x^) 



Proof of Theorem II. II 

We take ^ = 20e and a — 6^/2. Assume that L{h) is finite. By hypothesis and 
Theorem 4.1 xq is an approximate zero of fa- {ftXt)t&[a,b] denotes the unique 
lifting of the path ft corresponding to Co zero of fa associated to a;o. 
Induction hypothesis: 

£13/2 

-^tJ-{fti,x,)l3o{ft^,Xi) < a 

and furthermore, Xi is an approximate zero associated to Ca ■ 

The induction hypothesis holds by hypothesis at i = 0, so we assume it is verified 
up to step i. We are in the hypotheses of Lemma 5.1 for t = ti, t' — i^+i = 
min(s, s', h) and x — Xt- — xi. Thus, 

1)3/2 
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and iVjj {xi) is an approximate zero for /f._^i associated to Cti+i- 

Now we are in the hypotheses of Lemma 5.2 Thus, y = a^i+i picked at (1.3) 
satisfies the induction hypothesis. 

In order to bound the number of iterations, we remark that at each step z, one 
of the following alternatives is true: 

(1) This is the last step: i^+i — b. 

(2) Condition (1.1) is true. In that case, we are under the hypotheses of 
Lemmal5.3l 



(3) Condition (1.2 1 is true. Then we are under the hypotheses of Lemma 5.4 
Therefore, we may infer that at each non-terminal step. 



i+l) 



Liu) > 



13Z3'V2 ■ 



Therefore, there can be no more than 260i^ ^L{b)D^/^ non-terminal steps. There 
is only one terminal step, so the total number of steps is at most 



1 + 2mC'^L{b)D^''^ = 1 -I- Qme-'^L{b)D^/'^ 
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